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ABSTRACT 

We have used the ROSAT PSPC to study the properties of a sample of 24 X-ray 
bright galaxy groups, representing the largest sample examined in detail to date. Hot 
plasma models are fitted to the spectral data to derive temperatures, and modified 
King models are used to characterise the surface brightness profiles. 

In agreement with previous work, we find evidence for the presence of two com- 
ponents in the surface brightness profiles. The extended component is generally found 
to be much flatter than that observed in galaxy clusters, and there is evidence that 
the profiles follow a trend with system mass. We derive relationships between X-ray 
luminosity, temperature and optical velocity dispersion. The relation between X-ray 
luminosity and temperature is found to be Lx oc T 4 9 , which is significantly steeper 
than the same relation in galaxy clusters. These results are in good agreement with 
preheating models, in which galaxy winds raise the internal energy of the gas, inhibit- 
ing its collapse into the shallow potential wells of poor systems. 
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1 INTRODUCTION 

The majority of galaxies in the universe are found in galaxy 
groups (Tully 1987). These collections of between 3 and 
about 30 galaxies trace large scale structure (Ramella, Geller 
& Huchra 1990) and probably contain a large fraction of 
the total baryonic mass in the universe (Fukugita, Hogan & 
Peebles 1998). However despite their abundance and impor- 
tance, galaxy groups have received relatively little attention 
until recently. The main problem has been the identification 
of the groups themselves. Even when redshift information is 
available, it is difficult to identify whether a group is truly 
bound, due to the problems of small number statistics and 
chance superpositions. In contrast, galaxy clusters which are 
easier to identify due to the larger number of members, have 
been extensively studied. 

The detection of extended X-ray emission from hot gas 
in the group potential well provides the best evidence that 
a group is truly gravitationally bound. The study of this 
hot intragroup gas can provide important insights into the 
evolution and dynamics of the group and its member galax- 
ies. Samples of X-ray bright groups were originally studied 
using the Einstein satellite (e.g. Price et al. 1991), but the 
introduction of the ROSAT satellite with its improved sen- 
sitivity and resolution, allowed a more thorough analysis of 
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these systems. Since the ROSAT PSPC was first used to 
study X-ray bright groups (Mulchaey et al. 1993; Ponman 
& Bertram 1993) a number of collections of groups have been 
studied (e.g. Doe et al. 1995; Ponman et al. 1996; Burns et al. 
1996; Mulchaey et al. 1996; Mahdavi et al. 1997; Mulchaey 
& Zabludoff 1998). However, none of these studies provides 
a uniform, detailed analysis of a reasonable sized sample of 
groups, based on high quality data. The largest samples have 
all been based on ROSAT All Sky Survey (RASS) data, in 
which case properties other than the luminosities are diffi- 
cult to determine, due to poor statistics resulting from the 
short exposures. 

The study of Ponman et al. (1996) used a mixture 
of RASS and pointed data and identified 22 X-ray bright 
groups. These were all compact groups from the catalogue 
of Hickson (1982). Such compact groups have the advantage 
that they can be easily identified on the sky due to the high 
projected over-densities of galaxies within them, but may be 
unrepresentative of groups as a whole. The X-ray properties 
of these Hickson compact groups (HCGs) showed systematic 
departures from those of clusters, leading to the suggestion 
that they might be displaying the marks of energy injection 
into the intergalactic medium due to galaxy winds. 

Mulchaey & Zabludoff (1998) (henceforth MZ98) used 
pointed PSPC data to study groups of both loose and com- 
pact morphology. With a sample of only nine groups they 
were unable to derive reliable statistical results, however 
they found that properties such as the L : T relation and 
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surface brightness slope were indistinguishable from those 
of clusters, in contradiction to the results of Ponman et al. 
(1996). If this is true, it suggests a fundamental difference 
between the properties of loose and compact groups. The 
main aim of the present work is to establish the X-ray prop- 
erties of loose groups by means of a careful and uniform 
study of a larger sample of systems, and to establish whether 
they do, in fact, differ from compact groups in their X-ray 
properties. To allow direct comparison with the results of 
MZ98, their sample has been included within ours. 

In § 5 we describe the sample selection and initial identi- 
fication of the X-ray bright groups. The spectral and spatial 
analyses of the X-ray emission are described in § |^ and § ^ 
Results of the analysis, including correlations between the 
derived parameters, are presented in § |^. These results are 
compared with those of MZ98 in § |^, and discussed in § [?j 
Finally, our conclusions are summarized in § ^[ Throughout 
this paper we use Ho = 50 km s^ 1 Mpc -1 . 



2 SAMPLE SELECTION AND DATA 
REDUCTION 

The primary aim of this work is to study the properties of 
a number of X-ray bright groups, as such it was necessary 
to initially compile such a sample. Three different sources 
were used for this purpose, the optical catalogue of Nolthe- 
nius (1993), the sample of Ledlow et al. (1996) and the X- 
ray bright groups from MZ98. The catalogue of Nolthenius 
(1993) contains 173 groups, with three or more members, 
selected from the CfAl galaxy redshift catalogue using a 
friends of friends algorithm with a density enhancement of 
15. The Ledlow et al. (1996) sample contains 71 groups se- 
lected from the Zwicky Catalogue of Galaxies and Clusters 
of Galaxies, using a friends of friends algorithm with a sur- 
face density enhancement of 46.4, each group having at least 
4 members and galactic latitude |b| > 30°. 

Cross-correlation of the Nolthenius (1993) and Led- 
low et al. (1996) samples with the ROSAT observing log, 
identified groups which had been observed by the ROSAT 
PSPC during its programme of pointed observations. We 
further restricted ourselves to groups which had been ob- 
served within 20' of the centre of the PSPC. The nine X-ray 
bright groups from MZ98 were all known to have been ob- 
served by the ROSAT PSPC and were added to the sample. 
Groups identified as being part of known bright galaxy clus- 
ters such as Coma were also excluded at this stage. This 
resulted in a potential sample of 37 galaxy groups, which 
are listed in Table |l| 

Before the X-ray data can be used it is necessary to 
identify and exclude sources of contamination. These include 
particle events and solar X-ray emission scattered from the 
Earth's atmosphere into the telescope. Detectors on board 
the spacecraft identify and exclude over 99% of the parti- 
cle events that would be recorded as X-rays. These particle 
events are recorded as the master veto rate. At values of 
the master veto rate of above 170 count s _1 the contamina- 
tion by particles is significant, and these times are excluded 
from our analysis. Reflected solar X-rays can be identified 
by an increase in the total X-ray event rate. To remove this 
contamination, times where the total event rate deviated by 



more than 2a from the mean were excluded. Typically this 
resulted in the removal of a few percent of each observation. 

A standard reduction of the data was then carried out 
to produce an image and background for each group. The 
statistical significance of any emission within distances of 
50 kpc and 200 kpc from the optical centre of each of the 
groups was then calculated. This was used along with a 
smoothed image and a profile of the group to identify the 
presence of extended emission above a 5a detection thresh- 
old. It was also apparent that in a few cases diffuse X-ray 
emission was centred on a galaxy within the PSPC ring, 
even though the catalogued optical centre of the group was 
outside the ring. These groups were also included, and are 
identified with an asterisk in Table [l]. This resulted in a final 
sample of 24 X-ray bright galaxy groups, which are identified 
in Table [l]. As the table shows, those systems span a consid- 
erable range in catalogued optical richness (N ga i = 3 — 45). 
A more reliable measure of the total mass of each group is 
given by the X-ray temperatures derived below. Our sample 
should not be regarded as being statistically complete in any 
way, but rather a reasonably representative sample of X-ray 
bright groups. 



3 SPECTRAL ANALYSIS 

Events surviving the initial screening process were binned 
into a 3-dimensional x, y, Energy data cube. An estimate of 
the background was generated from an annulus at r=0.6- 
0.7°with the PSPC support spokes removed. The dataset 
was then background subtracted, and point sources identi- 
fied using a maximum likelihood source searching program. 
Point sources within the background annulus were removed 
to 1.2 times the 95% radius for 0.5 keV photons. The back- 
ground was then recalculated and the image once again 
searched for point sources. Other more extended sources, 
such as background galaxy clusters not associated with the 
group emission, were also manually identified and excluded 
at this point. This process of identifying and removing point 
sources to produce a better estimate of the background was 
repeated until the same number of point sources was iden- 
tified each time. Typically this took 4-5 iterations for each 
dataset. 

The final background subtracted data were then cor- 
rected for dead time effects and vignetting, and then divided 
by the effective exposure time to give a map of spectral flux. 
A circular region around each of the groups was used to ex- 
tract a spectrum. The size of this region was determined by 
examination of a smoothed image and a surface brightness 
profile of the group. The region was selected to include all 
the emission that could be observed in the smoothed image 
and profile; its size for each of the groups is shown in Ta- 
ble [l. Point sources, and other sources as identified above, 
were removed from the spectral image, along with the sup- 
port structure and the data outside the radius of interest. 
The spectrum for each group was then obtained by collaps- 
ing the spectral image along the x and y axes. 

Each spectrum was fitted with a MEKAL hot plasma 
model (Mewe, Lemen & van den Oord 1986) with a hydrogen 
absorbing column frozen at a value determined from radio 
surveys (Stark et al. 1992). For two of the groups it was 
also necessary to fix the abundance to obtain a sensible fit. 
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Table 1. Listed are the groups in which a search for extended X-ray emission was carried out. Groups with a 1 in the comments 
column have properties taken from Nolthcnius (1993), those identified with a 2 are from Ledlow et al. (1996) and those marked with 
a 3 are from Mulchaey & Zabludoff (1998). Asterisks indicate groups in which emission was identified within the PSPC support 
ring, but whose catalogued optical positions were outside the ring. Groups with detected X-ray emission are listed in the top half 
of the table along with the radius to which emission was observed. Groups that were not used are given in the lower region of the 
table along with the reason for exclusion. 
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A value of 0.3 solar was used for this purpose. In this way 
we derived temperature, abundance and bolometric flux for 
each group. 

For hot spectra, the limited spectral band of ROSAT 
makes temperature determination subject to systematic er- 
rors in the high energy response of the PSPC, and there 
is evidence that ROSAT temperatures are systematically 
lower than those from hard X-ray instruments such as Ginga 
and ASCA. A comparison of ROSAT and ASCA tempera- 
tures by Hwang et al. 1999 showed that this temperature 
bias amounts to ~ 30% in hot systems, but that there is 
no evidence of any systematic offset below an ASCA de- 
rived temperature of 2keV, where the ROSAT band covers 
the spectrum adequately. We therefore expect the ROSAT- 
determined temperatures for the systems in our sample 
(which have T < 1.7 keV) to be free from serious bias. 

The distribution of group temperatures in the present 
sample occupies a rather small range around 1 keV. For 



the 22 groups in which metallicities were derived, the over- 
all weighted mean metallicity is 0.19 ± 0.01 solar, whilst 
the median is 0.42 solar. A trend is observed in clusters for 
higher metallicity in lower temperature systems (Arnaud 
1994). This would lead one to expect a typical metallicity 
of ~ 0.6 solar in systems with T ~ 1 keV. However there 
is evidence that there may be abundance gradients in cool 
clusters which result in an increased abundance at the cen- 
tre (e.g. Xu et al. 1997; Ikebe et al. 1997; Fukazawa et al. 
1998; Finoguenov & Ponman 1999) and could account for 
the observed trend. In any case, results obtained here for 
the group metallicities must be viewed with caution, since 
ROSAT is unable to resolve individual emission lines, and 
metallicities can be strongly biased when a variable temper- 
ature plasma is fitted with an isothermal model (e.g. Buote 
& Fabian 1998, Finoguenov & Ponman 1999). 

For each of the groups in the sample we also derived 
simple projected temperature profiles. In each case spectra 



© 1998 RAS, MNRAS 000, 



4 Stephen F. Helsdon and Trevor J. Ponman 



in several annuli were extracted and fitted with a MEKAL 
model as described above. In each annulus the hydrogen col- 
umn and abundance were frozen at the global values. The re- 
sulting temperature profiles of all groups are shown in Fig.jlJ 
Some of the profiles shown are not particularly informative 
due to a combination of large errors on the temperatures 
and a small number of annuli. However it is clear that ap- 
proximately half of the groups show evidence of a tempera- 
ture drop in the central regions, indicating the presence of 
a cooler component. Also, approximately half of the profiles 
show evidence of a decline in temperature at large radii. 



4 SURFACE BRIGHTNESS PROFILES AND 
GROUP LUMINOSITIES 

Observations of galaxy clusters across a wide range in virial 
temperature appear to indicate a flattening of the profiles 
in lower mass systems (Arnaud & Evrard 1999; Ponman, 
Cannon & Navarro 1999) - a result consistent with expecta- 
tions if the intergalactic gas has been subject to preheating 
by galaxy winds (Metzler & Evrard 1999; Cavaliere et al. 
1999). However, MZ98 found for their sample of groups that 
surface brightness profiles did not differ significantly from 
those in clusters, once the presence of central components 
was properly allowed for. We set out to examine the sur- 
face brightness profiles of our group sample in an attempt 
to resolve this issue. 

Following initial reduction, an image was extracted in 
the 0.5 - 2 keV band, and corrected for vignetting using an 
energy dependent exposure map (see Snowden et al. (1994) 
for description) . Point sources identified in the spectral anal- 
ysis were removed from the image along with any other un- 
related extended sources. Only the data within the region 
from which each group spectrum was extracted, were used 
for the spatial analysis. It has been shown that the centroid 
of the X-ray emission often lies at the position of the bright- 
est group galaxy (MZ98), and as such any emission centred 
on this galaxy may be associated with the group potential as 
a whole. For this reason any source associated with the cen- 
tre of the X-ray emission was not removed. Use of the energy 
dependent exposure map to correct for vignetting, results in 
a constant background level across the image, therefore a 
flat background was also determined and subtracted from 
the data. 

For each group the 2-dimensional surface brightness 
profile was modelled with a modified King function (or '/3- 
profile')of the form: 

5( r )=So(l + (r/w) 2 )- 3 ^+ a5 

Models were convolved with the PSPC point spread 
function at an energy determined from the mean photon 
energy of the group spectrum, and fitted to the data. The 
free parameters were the central surface brightness So, the 
core radius r core , the index f3ju and the x and y position of 
the centre of the emission. Both spherical and elliptical fits 
were carried out on the data, with the major to minor axis 
ratio and the position angle being extra free parameters in 
the elliptical fits. 

The use of 2-dimensional datasets to fit the surface 
brightness distribution results in a low number of counts in 



many of the data bins. Under these conditions chi-squared 
(x 2 ) fitting performs poorly, so maximum likelihood fitting, 
using the Cash statistic, was used instead. The Cash statis- 
tic (Cash 1979) is defined as — 21nL where L is the likelihood 
function (in this case derived from the Poisson distribution). 
Thus the most likely model has a minimum Cash statistic. 
Differences in the Cash statistic are \ 2 distributed, so con- 
fidence intervals may be calculated in the same way as for a 
conventional \ 2 fit. 

Unfortunately the Cash statistic by itself gives no indi- 
cation of the quality of a fit; hence it was necessary to obtain 
some other estimate of the fit quality. A Monte Carlo ap- 
proach was used, in which the best fit model was used to 
simulate 1000 images of the group. Poisson noise was added 
to each of these images, and they were then compared to the 
original model, and the Cash statistic for each image deter- 
mined. Thus, for a particular model we were able to obtain a 
distribution showing the range of Cash values expected for 
datasets generated from this model. A Gaussian was then 
fitted to this distribution to obtain the width and central 
value. By comparing the Cash statistic for the real dataset 
with this distribution, it was possible to determine the prob- 
ability that the model could have produced the data. This 
probability is recorded in Table ^[ as the number of standard 
deviations that the real value lies from the centre point of 
the distribution. If the value of the real Cash statistic lay 
more than 2a from the peak of the distribution then the fit 
was regarded as 'poor'. 

As can be seen in Table ^| the single-component fits pro- 
vide an adequate description of the data in a few cases. How- 
ever for most groups the single-component fits are poor. It 
has been suggested that there are typically two components 
in the surface brightness profiles of galaxy groups (MZ98) , a 
central component associated with a central galaxy, cooling 
flow or AGN, and a more extended component associated 
with the group potential. To check this, models compris- 
ing of two superposed /9-profiles were also fitted to those 
datasets with poor single-component fits and greater than 
~ 900 total counts. Below this number of source counts, 
statistics were found to be too poor to constrain the more 
complicated two-component models. To limit the number of 
free parameters, the central component was constrained to 
be spherical while the outer component was allowed to vary 
in ellipticity. 

In three of the groups (NGC4065, NGC4073 and 
NGC7619) the emission was bimodal, so that the two- 
component models fitted with the centres of the two compo- 
nents significantly offset from one another (e.g. see FigM. 
As a result, it is not sensible to define one component as 
extended, and the other as the central component. In these 
cases both of the components were constrained to be spher- 
ical. 

The fitted parameters of the two-component King pro- 
files are also shown in Table ^, along with an estimate of 
the goodness of fit. The errors quoted are la for one inter- 
esting parameter. Note that these errors are only reliable 
for reasonable fits (see final column in Table |2|). The best 
fitting surface brightness profiles were also used to correct 
the derived group fluxes for the diffuse emission lost when 
point sources are removed. A model image for the group was 
produced, and from the ratio of the number of counts in the 
model image to that in the same image with 'holes' punched 
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Figure 1. Temperature profiles from an annular analysis of the 24 groups in this sample. A number of the groups show a temperature 
drop in the central regions and/or at large radii. 

© 1998 RAS, MNRAS 000, 



6 Stephen F. Helsdon and Trevor J. Ponman 



Table 2. Results of the surface brightness fits for elliptical and two-component models. If two models are shown for a group, the first is 
the elliptical model and the second the the two-component model. Models marked with an * are groups in which two separate centres of 
emission could be observed. The goodness of fit is as described in the main body of the text. All errors are lcr for one interesting parameter. 



Extended component Central component 



Group 


Pa 


t 


Core radius 
(arcmin) 


Axis ratio 


Position angle 
(degrees) 


Core radius 


Pfit 


Goodness 
oi nt. 


IN Li V_/ 


OlO 


1 Q7 
l.O f 


± 


0.36 


0.42 ± 0.10 


l.Uo 


it U.lo 


141 ± 65 


_ 


_ 


-0.5 


Mpp 
IN KjKj 


OOO 


n ^9 
u.ooz 


Zt 


0.003 


0.43 ± 0.06 




□z u.uo 


166 ± 4 


_ 


_ 


8.2 






U.4o 


Zt 


0.02 


6.9 ± 0.9 


1 1 o 

i. ty 


4- n r\A 

± U.U4 


156 ± 6 


0.01 ± 0.03 


0.48 zt 0.02 


5.2 


In LtO 


oz4 


U.40 


1 

Zt 


0.01 


0.01 ± 0.01 


1 A 1 
1.41 


± U.ZZ 


143 ± 16 


_ 


_ 


18.2 


IN 


K Q "3 


U.40Z 


Zt 


0.005 


0.40 ± 0.03 


1 fin 
l.DU 


± U.UO 


77 ± 3 


_ 


_ 


14.2 






U. / O 


1 

Zt 


0.14 


10.2 ± 2.8 


1 CO 


I n r\A 
± U.U4 


49 ± 2 


0.02 ± 0.02 


0.54 zt 0.02 


-0.23 


Mpp 
IN Kj v_ 


1 1 "1 
Ml 


U.400 


1 

Zt 


0.008 


0.37 ± 0.05 


1 Qc; 

i-.OO 


-U n no 
it u.uy 


74 ± 6 


_ 


_ 


11.2 






U.oy 1 


_1_ 

Zt 


0.009 


0.10 ± 0.13 


L.Z<± 


i n 11 

it U.ll 


174 ± 11 


0.17 zt 0.09 


0.9 zt 0.3 


0.20 


M p p 

IN 


lOo / 


U.4f 


_1_ 

Zt 


0.06 


0.34 ± 0.27 


1 .1 


± u.o 


35 ± 22 


_ 


_ 


-0.1 


Mpp 
IN 


ZOOo 


u.ooy 


_l_ 

Zt 


0.003 


0.01 ± 0.02 


l.ZO 


It u.u^t 


88 ± 7 






1.6 






n /inn 

U.4UU 


1 

Zt 


0.004 


2.6 ± 0.6 


l.ol 


± U.UD 


95 ± 7 


0.2 Fixed 


1.0 Fixed 


1.54 


Mpp 
IN 


ouyi 


U.OU 


_1_ 

zt 


0.02 


0.48 ± 0.06 




4- HQ 

it u.uy 


14 ± 4 


_ 


_ 


10.0 






U.41 


zt 


0.02 


0.1 ± 0.05 


l.O 


it l.O 


184 ± 4 


0.3 ± 0.07 


0.66 zt 0.05 


0.82 


MPP 
IN Li V_/ 


oOU / 


n 

yj.oz 


_1_ 


0.18 


4.8 ± 1.2 


Q 1 
O.'l 


4- n v 


20 ± 3 


_ 


_ 


3.2 








ZL 


0.04 


0.28 ± 0.16 


i n 

1 A) 


Fixed 


0.0 Fixed 


0.01 ± 0.04 


0.38 zt 0.03 


0.63 


MPP 
IN y i v_ 


oODO 


U.4y 


1 

zt 


0.03 


0.13 ± 0.12 


"1 A 
1 . 1 


± (J.Z 


63 ± 19 


_ 


_ 


0.2 


wpp 

IN 




U.I / 


_L 
ZL 


0.04 


4.1 ± 0.6 


. 


4- n a 
It u.^t 


7 ± 1 


_ 


_ 


8.9 






U.441 


_L 
ZL 


0.01 


0.05 ± 0.06 


i n 

l.U 


Fixed 


0.0 Fixed 


4.3 zt 1.4 


0.8 zt 0.2 


2.49 


MPP 

IN y i v_ 


4U ( o 


U.4ol 


1 

Zt 


0.002 


0.10 ± 0.01 


L.ZO 


4- n no 

± U.UZ 


103 ± 3 


_ 


_ 


2.5 






U.10 


_L 
ZL 


0.01 


2.14 ± 0.26 


1 9S 
l.ZO 


4- n c\a 
It u.u^t 


102 ± 4 


0.34 ± 0.05 


0.73 zt 0.06 


1.64 


MPP 
IN Li V_/ 


/I 9K1 

4zoi 


U.440 


1 

Zt 


0.004 


0.01 ± 0.01 


1 *3 A 
1.04 


4- n no 
it U.uy 


47 ± 5 


_ 


_ 


13.6 






U.oO 


_L 
Zt 


0.03 


3.3 ± 1.0 


1 99 
l.ZZ 


J. fl 19 
It U.lo 


156 ± 17 


0.28 zt 0.01 


1.0 Fixed 


-1.20 


MPP 
IN y i v_ 


/I Q9^ 


n An 

U.DU 


1 

Zt 


0.01 


0.32 ± 0.03 


1 17 
1.1/ 


it U.UO 


11 ± 8 


_ 


_ 


1.3 


MPP 
IN Kj v_ 


40oO 


U.4 < 


_l_ 

Zt 


0.003 


0.268 ± 0.014 


i no 
i.uy 


4- n fie 

it U.Uo 


17 ± 7 






17.0 






0.373 


± 


0.008 


0.013 ± 0.006 


1.15 


ib 0.06 


167 ± 9 


0.92 zt 0.11 


0.81 zt 0.06 


6.39 


NGC 


4761 


0.502 


± 


0.005 


0.30 ± 0.02 


1.33 


ib 0.04 


16 ± 3 


_ 


_ 


19.3 






0.364 


zt 


0.006 


0.10 ± 0.04 


1.17 


ib 0.09 


72 ± 14 


0.49 zt 0.06 


0.85 zt 0.04 


1.28 


NGC 


5129 


0.44 


± 


0.01 


0.15 ± 0.06 


1.30 


ib 0.16 


28 ± 13 






-0.7 


NGC 


5171 


0.34 


± 


0.03 


0.84 ± 0.71 


2.7 


ib 0.9 


171 ± 6 






0.8 


NGC 


5353 


0.58 


± 


0.03 


1.35 ± 0.17 


1.34 


ib 0.11 


129 ± 7 






71.0 






0.44 


± 


0.02 


0.29 ± 0.13 


1.8 


ib 0.2 


25 ± 5 


0.01 zt 0.01 


1.0 Fixed 


6.39 


NGC 


5846 


0.66 


± 


0.02 


1.43 ± 0.08 


1.70 


ib 0.05 


246 ± 2 






106.9 






0.58 


± 


0.01 


0.84 ± 0.07 


1.16 


ib 0.03 


47 ± 7 


0.27 zt 0.04 


1.0 Fixed 


3.48 


NGC 


6338 


0.423 


± 


0.004 


0.06 ± 0.01 


1.24 


ib 0.06 


24 ± 7 






21.9 






0.52 


± 


0.04 


2.6 ± 0.4 


1.29 


ib 0.07 


23 ± 6 


0.42 zt 0.03 


1.0 Fixed 


1.40 


NGC 


7176 


1.07 


± 


0.29 


8.0 ± 4.2 


1.6 


ib 0.2 


139 ± 10 






2.3 


NGC 


7619 


0.458 


± 


0.006 


1.42 ± 0.08 


1.93 


ib 0.07 


210 zfc 1 






22.7 


* 




0.78 


± 


0.08 


0.26 ± 0.05 


1.0 


Fixed 


0.0 Fixed 


0.01 zt 0.01 


0.40 zt 0.01 


6.27 


NGC 


7777 


0.35 


± 


0.02 


0.01 ± 0.01 


1.0 


Fixed 


0.0 Fixed 






1.1 



at the positions of detected sources, a correction factor for 
the fluxes was obtained. The luminosity of each group was 
then calculated using distances corrected for infall to Virgo 
and the Great Attractor (Fixsen et al. 1996; Burstein 1990), 
which are listed in Table 

In the case of groups with a detected central component, 
we checked for the possibility that this might arise from a nu- 
clear source in the central galaxy. Fits with Gaussian models 
for the central component show that it is extended at > 99% 
confidence in all cases except NGC 5353, where statistics are 
too poor to constrain the extent of the central source. For 
each of these systems a search for radio sources associated 
with the brightest group galaxies was carried out using NED 
and the Burns et al. 1987 radio survey of groups, which has 
some overlap with this sample. This search identified radio 
sources associated with six of the brightest group galaxies: 
NGC 383, NGC 741, NGC 4261, NGC 4636, NGC 5353 and 



NGC 6338. Hwang et al. 1999 have studied three of these us- 
ing ASCA spectra. Two showed no significant improvement 
in fit statistic when a powerlaw component was added to 
the spectral model. For the third, NGC 6338, Hwang et al. 
1999 find evidence that there may be contamination by an 
AGN, although in our data the spatial extent of the central 
component in this, and almost all the other systems rules 
out a large AGN component. The only system that may be 
contaminated (as indicated by the spatial extent of the cen- 
tral component) is NGC 5353. For this system we fitted the 
spectral data for this group with an added power-law com- 
ponent of index 1.7. We then calculated the relative contri- 
butions from the power-law and hot plasma components in 
the ROSAT band. This showed that even for the 90% upper 
limit of the power-law component the emission was domi- 
nated by the hot plasma component. Our conclusion from 
these spatial and spectral studies is that any AGN contri- 
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Figure 2. Contours of adaptively smoothed X-ray emission from 
the group NGC4065, overlaid on an optical image. It is clear that 
there are two distinct centres of emission in this system. 



bution to the central components in these systems appears 
to be minor. 

We were also interested in the way in which the lu- 
minosity of groups varies with radius. The model images 
were therefore used to calculate luminosities within radii of 
200 kpc, 500 kpc, 1/3 of the virial radius and the virial 
radius (Rv)- Note that Rv lies well beyond the radius to 
which significant X-ray emission can actually be detected in 
our data, in almost all cases. It can be seen in Table ^ that 
the two-component models provide good descriptions of the 
data in the majority of cases. However even in the cases 
where the two-component fit is not acceptable it is signifi- 
cantly better than the single-component model, thus where 
possible, the two-component models are used for the pur- 
poses of calculating the effects of extrapolating to different 
radii. The virial radii of the groups were determined using 
the relation obtained from simulations by Navarro, Frenk & 
White (1995). This is given (for H = 50 km s" 1 Mpc -1 ) 
by, 

^ = 2.57(33^)4 Mpc. 

Luminosities and temperatures derived in this study are 
generally similar to those from earlier studies (Doe et al. 
1995; Ponman et al. 1996; Mulchaey et al. 1996; Burns et al. 
1996; Mulchaey & Zabludoff 1998). The small differences in 
the luminosities of groups common to both this sample and 
that of Burns et al. (1996) are most likely primarily due to 
the fact that Burns et al. (1996) use a spectral model with 
a temperature of 1 keV to derive all luminosities, whereas 
the luminosities derived here use fitted spectral models, and 
thus should be more reliable. 



NGC 2563 




Radius (arcmin) 



NGC 3091 




Radius (arcmin) 



Figure 3. ID surface brightness profiles for the groups NGC 2563 
and NGC 3091. The overall best fit two component models are 
shown as the solid lines, with the dotted line representing the 
central component and the dashed line the extended. Data points 
are shown as crosses. For comparison the single component ellip- 
tical models are marked as the dot-dashed lines. This is steeper 
than the extended component in NGC 3091 but slightly flatter 
than it in NGC 2563. r core marks the core radius of the extended 
component and r x the crossover radius, as defined in the text. 

5 RESULTS 

Throughout the following sections, the luminosities quoted 
are extracted from within the radius given in Table pi Cor- 
rections for removed point sources have been made using the 
best model derived for each group; either two-component or 
(elliptical) single-component. 

5.1 X-ray profiles 

The surface brightness profiles for our 24 systems break 
down into 12 two-component, 9 single-component and 3 bi- 
modal cases. However, note that the nine single-component 
systems include the eight groups with the lowest source 
counts in the sample, so it is likely that the majority of 
these single-component fits appear to be adequate only be- 
cause of poor statistics. Two examples of radial profiles are 
shown in Fig.^[ These ID profiles only give an approximate 
representation of our 2D models, but the centres of the two 
components almost coincide in the two cases shown, and 
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Table 3. Results of the spectral fitting are shown along with the derived luminosities for each of the groups. The distance to each of 
the groups is calculated after allowances for infall to Virgo and the Great Attractor. The luminosities are those derived within the radius 
to which emission could be observed. The final column shows what fraction of the luminosity as extrapolated to the virial luminosity is 
observed within the stated radius. All errors are la. 
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* NGC 315 fits with a high Pfu and it is possible that for this group the emission may be due to a extensive elliptical galaxy halo rather 

than genuine group emission. 



profiles for both data and model components have been de- 
rived about the centre of the more compact component. A 
distinct 'shoulder' in the observed profile indicates the need 
for two components in the model, as noted by MZ98. 

The median value of f3[u obtained for the extended 
group component in our sample (from two-component fits 
where available, or else single-component fits), is f3f it = 0.46, 
and the weighted mean is 0.42 ± 0.06 (where the error is de- 
rived from the scatter of the values about the mean). This 
value of Pfu can be compared with the typical value for rich 
clusters, (3 fit « 2/3 (Arnaud & Evrard 1999; Mohr, Math- 
iesen & Evrard 1999)), indicating that the surface brightness 
profiles of the groups in this sample are generally signifi- 
cantly flatter than those of clusters. 

In the case of a number of the groups, the core radius 
derived for the extended group component is smaller than 
the resolution of the ROSAT PSPC. Such a small core ra- 
dius means essentially that the group emission has been fit- 
ted with a power law model. In such cases, the derived core 
radii are unreliable, particularly if an additional central com- 
ponent is present. To investigate the effect of a small core 
radius on the other derived parameters, in particular the 
slope, Pfu, we varied the core radii in a number of groups 
which fitted with two components and a small core radius 
for the extended component. It was found that varying the 
core radii between 0.1 and 1.0 arcmin typically changes the 
value of Pfa by less than 5%. Values of the index and core 
radius of the central component also varied only within lcr 



of their best fit values. Hence uncertainties in core radius in 
such cases do not seriously compromise our results for other 
parameters. 

We define the 'cross-over radius', r x , to be the radius 
within which the central component dominates the surface 
brightness. We derived values for r x for the twelve groups 
for which two-component fits were available, but the emis- 
sion was not bimodal. The mean cross-over radius for these 
systems was r x = 35 ± 6 kpc. Groups with two-component 
profiles which have r core < r x are deemed to have poorly 
determined core radii. In order to gain some insight into the 
typical core radii of galaxy groups, the median value was de- 
termined, using r x as an upper limit for those groups with 
r CO re < r x . Under these assumptions the median core radius 
of the twelve groups was found to be 60 kpc. 

The relationship between the integrated temperature 
of the intragroup gas and the best obtained value 
for each of the groups is plotted in Fig^. Also shown are 
cluster data from Arnaud & Evrard (1999) (data points 
with circles). The group data are split into two categories: 
single-component (plain crosses) and the extended compo- 
nent from the two-component fits (points with square in 
centre). As can be seen, the general trend in clusters is for 
Pfit to drop with decreasing temperature. In the region of 
the graph containing the group data it is clear that the ma- 
jority of the groups have low (3f it values, but there is also 
a large amount of scatter, in particular amongst the groups 
with single-component fits. 
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Figure 4. The relationship between f3fu and temperature for the whole group sample compared to cluster data from Arnaud & Evrard 
(1999). Values from groups with single-component fits are shown as plain crosses, those from two-component fits as crosses with central 
squares. The cluster data are marked as crosses with central circles. 



Fig.y also includes the three bimodal groups. Exclud- 
ing these three groups, the single-component fits and the 
two-component fits whose quality of fit (from the Monte 
Carlo simulations) is poor, greatly reduces the scatter in 
the group results. The outcome is a much clearer trend in 
(3fit with temperature, as can be seen in Fig.[]. The com- 
bined groups and cluster data are significantly correlated 
with a Kendall's rank correlation coefficient (a distribution 
free test for correlation) of K=4.05 (P=0. 00006 of chance 
occurrence). The one group point that conflicts with the 
general trend (NGC533) is, in fact, the only group in the 
sample which has a flatter j3fu value for its central com- 
ponent than for the extended component. This means that 
the central component has a significant effect beyond the 
central region. Hence the shape of this component could af- 
fect the parameters obtained for the extended component. 
To test this, we refitted the surface brightness profile with 
a Gaussian model for the central component, in place of the 
previous King model. The j3fu value for the extended com- 
ponent changed markedly, and the new value is denoted by 
the triangle in Fig.^. As can be seen, this point is now much 
closer to the trend described by the other groups. 



5.2 Luminosity, temperature and velocity 
dispersion 

5.2.1 X-ray luminosities 

Bolornetric luminosities for each group, derived from within 
the extraction radius as described in section 3, are given in 
Table ^[ along with best fit spectral properties. The tabu- 
lated luminosities are those of the intragroup gas only. Er- 
rors on the luminosities are derived from Poisson errors on 



the data combined with errors arising from uncertainties in 
the fitted surface brightness profiles, which are used to cor- 
rect for flux lost where contaminating sources have been 
excluded. 

The flat surface brightness profiles of groups imply that 
a significant fraction of their luminosity derives from large 
radii. To quantify this, we used our best fit surface bright- 
ness models to derive bolornetric luminosities extrapolated 
to Rv, and the fraction of this luminosity represented by 
the luminosity derived from within the extraction radius is 
shown for each group in Table ^| This may be as low as 
~30% in some cases. 

The effect of scaling the luminosities to different radii 
is shown in more detail in Fig.^J. This analysis is based on 
the eight systems with well-fitting two-component profiles. 
These have been binned into three temperature bins to re- 
duce fluctuations from system to system and show trends 
more clearly. The luminosity as a fraction of that within Rv 
is shown at three radii for the systems within each temper- 
ature bin. Points marked by triangles (dash-dot-dot line) 
show these ratios at a radius of 200 kpc, squares (dashed 
line) at the radius out to which emission could be detected, 
and crossed circles (solid line) at one third of the virial ra- 
dius of the group. As can be seen the luminosity is signifi- 
cantly underestimated in all cases. In particular, for groups 
measured to a fixed radius of 200 kpc, and for the coolest 
groups at the extraction radius, one may underestimate Lx 
by factor of more than two. 

5.2.2 Correlations 

The well-known relation between X-ray luminosity and tem- 
perature is apparent in our sample. The two parameters are 
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Figure 5. The relation between fifit and temperature for a subsample of the groups with the best two-component models, and the 
cluster data. Symbols are the same as in the previous figure. The triangle marks the new value of the one discrepant group value when 
its central component was refitted with a Gaussian. 
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Figure 6. Fraction of the total luminosity observed within three 
different radii, as deduced from best fitting surface brightness 
models, for systems of different temperature. Squares show the 
luminosity within the radius to which emission could be observed, 
circles the luminosity within Ry /3, and triangles the luminosity 
within a fixed radius of 200 kpc. 



significantly correlated (K=4.81, P<0. 00001) and the rela- 
tion between them is shown in FigQ Neither the errors on 
Lx or T are negligible, and a doubly weighted technique 
made available through the ODRPACK package was used in 



this and following plots to determine the best fit line, 
log L x = (42.98 ± 0.08) + (4.9 ± 0.8) log T . 

This relationship is marked with its la error bounds in 
Fig.[?j. A best fit to the cluster L : T relation has been derived 
by White, Jones & Forman (1997). They obtain logLx = 
42.67 + 2.98 log T, which is marked as the heavy dashed line 
in Fig. ^. This line is much flatter than the best trend fit for 
the loose groups. 

The luminosities used in this plot are those within the 
radius of extraction. Figj^ shows that at this radius the 
luminosities will be underestimated, with the effect being 
greatest in the smaller mass systems. This means that if 
luminosities extrapolated to the virial radius were used, 
the L : T slope should be slightly flatter. This is in- 
deed found to be the case, with a best fitting relation of 
logLx = (43. 17 ± 0.07) + (4.2 ±0.7) log T, although the dif- 
ference in slope from the previous relation is not formally 
significant. 

If galaxy systems scaled with mass in a self-similar way, 
then one would expect Lx oc T 2 . The cluster relation is 
steeper than this, and our result for groups is steeper still. 
However, the relationships derived by White et al. (1997) 
do not take into account the effects of cluster cooling flows, 
and recent work suggests that the L : T relation may be flat- 
tened towards Lx oc T 2 when the effects of cooling flows are 
allowed for (Allen & Fabian 1998; Markevitch 1998). Such a 
flattening of the relation for clusters would raise its extrap- 
olation at low temperatures, accentuating the disagreement 
with the low luminosities observed in groups. 

In Fig.[| velocity dispersion is plotted against the X-ray 
luminosity for our sample. A strong correlation can be seen 
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Figure 7. The relation between X-ray luminosity and temperature for our group sample. The solid line shows the best fit relation to 
our data, with one sigma error bounds marked by dotted lines. The extrapolation of the best fitting cluster relation (White et al. 1997) 
is shown as the dashed line. 
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Figure 8. The relationship between X-ray luminosity and group velocity dispersion, a. The best fit to the data is shown as the solid 
line with the one sigma error bounds marked by the dotted lines. The extrapolation of the cluster relation (White et al. 1997) is shown 
as the dashed line. 



between these two parameters (K=3.97, P=0. 00006). A re- 
gression line fitted to the data gives 

log L x = (31.3 ± 2.8) + (4.5 ±1.1) log a , 



which is marked in Fig.|8| with its la error bounds. This re- 
lationship is somewhat flatter than the cluster trend given 
by White et al. (1997) of logLx = 25.84 + 6.38 log a (bold 
dashed line in Figj^). Dell'Antonio, Geller & Fabricant 
(1994) found evidence that the L : a relation may flatten 
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below a ~ 300 km s _1 . However they did not remove the 
galaxy contribution from the X-ray emission, and suggest 
that their flattening may arise from the galaxy contribution 
becoming significant at low luminosities. This flattening has 
also been confirmed by Mahdavi et al. (f997). In the work 
presented here, contaminating sources were removed, but a 
flatter relation than clusters is still seen. Our result is actu- 
ally consistent with that expected from self similar scaling 
of clusters, i.e. Lx oc a 4 . However, errors are large and there 
is a good deal of scatter, so that the disagreement with the 
cluster result is not highly significant, and requires further 
confirmation. 

A strong correlation between a and T is shown in Figj^ 
(K=3.82, P=0.0001). A regression line fitted to the data 
gives 

logo- = (2.57 ±0.03) + (1.1 ±0.2) logT , 

which is shown in Fig|] with its la error bounds. Also shown 
in Fig|] is the line f3 spec = 1, where fl spec is defined as the 
ratio of the specific energy in the galaxies to that in the gas. 
As can be seen, this /3 spec = 1 line is flatter than the relation 
for the loose group sample. However it is interesting to note 
that the higher temperature groups appear to be consistent 
with /3 S pec = 1, while the lower temperature groups appear 
to fall well below this relation. The extension of the best fit 
relation for galaxy clusters as determined by White et al. 
(1997) is shown as the dashed line in Fig.^| This line, given 
by logo - = 2.53 + 0.6 log T, is also significantly flatter than 
the relation determined for the loose group sample. 

The unweighted mean value of f3 spec for our sample is 
0.86 ± 0.13. However, with one exception, it is clear that 
f3 3pec is decreasing in the lower temperature (i.e. lower mass) 
systems. These results are in good agreement with those 
of Bird, Mushotzky & Metzler (1995), who predict a trend 
towards lower f3 spsc in smaller systems. 

The one low temperature point (NGC 3607) that has a 
high velocity dispersion is also deviant in the L : a plot. Ex- 
amination of the group members reveals that, of the three 
catalogued members, one is a large angular distance from 
the remaining two, and has a large difference in recession 
velocity. Also there is a further bright galaxy at the redshift 
of the group, which is very close to two of the catalogued 
members. The recession velocity of this galaxy is between 
those of the two catalogued galaxies, and is almost certainly 
a group member, although it was not classified as such by 
Nolthenius (1993). These two effects combined indicate that 
the true velocity dispersion of the group is probably consid- 
erably lower than our estimate, which is taken from Nolthe- 
nius (1993). 



6 COMPARISON WITH MULCH AEY AND 
ZABLUDOFF 

As discussed in the introduction, we have included the X- 
ray bright systems studied by MZ98, in order to allow a 
direct comparison of our results with theirs. This is impor- 
tant, since our conclusions about /3fu, /3 sp ec and the L : T 
relation all differ from MZ98. In Table ^ we show the best 
fit parameters as determined by MZ98 for the groups that 
both they and we fit with two-component models (Note 



that they also fit two component models to NGC 4325 and 
NGC 5129, whilst we find that single component elliptical 
models provide an adequate representation of our data for 
these systems). Whilst we confirm their conclusion that two- 
component fits are required to adequately represent most 
systems, it can be seen there are some significant differences 
between the two sets of results. 

The fitting techniques used by MZ98 differ from those 
used in this work. Since they work with radial profiles, their 
fits are necessarily ID models, with both components cen- 
tred at the same point. Their method firstly involved exclud- 
ing the central region and fitting for the outer component 
only. The central component was then fitted with the ex- 
tended component fixed at the values derived from the pre- 
vious fit. Thus at no stage were the two components allowed 
to fit simultaneously. The 2D models fitted in this work al- 
low the positions of the two components to vary and also 
permit elliptical models to be used. Parameters for the two 
components were also optimised simultaneously. The lower 
number of counts in each bin forced us to use maximum 
likelihood fitting rather than x 2 fitting, but the quality of 
the fits were checked using the Monte Carlo approach as 
described above. 

To demonstrate the dangers of a ID approach to fit- 
ting the surface brightness profiles we simulated an image 
of a group, in which the outer component was elliptical (axis 
ratio=1.5), and offset a short distance (3 arcmin) from the 
central component. These values were chosen to construct a 
fairly elongated and offset system to make any biases more 
obvious. A 2D fit successfully recovered the slope of the outer 
component (/3^ t =0.4). We then attempted to fit the data us- 
ing a ID approach. We initially extracted a profile centered 
on the brightest point in the group (the central component). 
This gave a profile with a shoulder and a clear central ex- 
cess. This profile was fitted using QDP with a /3-profile plus 
a constant background. Initially we fitted to the full profile, 
giving a value of j3fa » 0.7. We then progressively excluded 
the central regions and refitted the data. The fitted value 
of Pfu rose to a peak of ~ 0.9 before dropping as a larger 
central region was excluded. Thus it is possible, with the ID 
approach used by MZ98, to significantly overestimate the 
true value of f5fu- 

To decide whether the models of MZ98 referred to in 
Table ^| provide an acceptable fit to our data, we carried out 
a series of two-component fits with the index and core radius 
frozen at the MZ98 values. The components were also con- 
strained to be circular and centred in the same place. The 
Cash values for these models were then compared to the 
best fitting values derived earlier. The differences between 
the Cash statistic values are shown in the final column of 
Table ^[ As can be seen, the models using the MZ98 param- 
eters generally fall well outside the 99% confidence regions 
of our best fitting models (which corresponds to AC=-20.1). 
Hence it appears that our more sophisticated models do rep- 
resent the data significantly better. 

The most important difference in the surface brightness 
results is apparent in the j3fu value of the extended compo- 
nent. MZ98 obtain values consistent with f3fu ~ 1 whereas 
the values obtained here mostly lie in the region 0.4-0.5, with 
a median value for the extended component of j3fu = 0.46. 

MZ98 obtain lower values of /3f it when fitting single- 
component models, but find that the extended components 
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Figure 9. The relationship between group velocity dispersion, <r and temperature. The best fit to the data is shown as the solid line 
with the one sigma error bounds marked by the dotted lines. The extrapolation of the cluster relation (White et al. 1997) is shown as 
the dashed line. The dot-dashed line shows the locus along which (3 spec = 1. 



Table 4. A comparison of the two-component models fitted by Mulchacy & Zabludoff (1998) with 
those from this work. (3 fa and core radius values of the extended component for both sets of models 
are listed. The final column gives the difference in Cash statistic between the two models as fitted 
to our data; the negative sign indicating that the model fitted here gives the better fit. 





M&Z 


M&Z 


This work 


This work 




Group 


Pfit 


core radius 


Pfit 


core radius 


AC 






(arcmin) 




(arcmin) 




NGC 533 


0.83 


8.15 


0.74 


10.2 


-122.5 


NGC 741 


1.00 


14.08 


0.39 


0.1 


-45.1 


NGC 2563 


0.86 


11.15 


0.40 


2.6 


-26.4 


NGC 3091 


0.68 


3.61 


0.41 


0.1 


-41.0 


NGC 4761 


0.63 


9.00 


0.36 


0.1 


-129.2 


NGC 5846 


0.83 


13.93 


0.58 


0.84 


-263.7 



fit with systematically higher (5f it when a second compo- 
nent is included (this sort of effect was reproduced in our 
simulations mentioned earlier). The same effect is noted for 
a sample of clusters by Mohr et al. (1999), who give a use- 
ful discussion of the effect. Since core radius and f3fu are 
strongly positively correlated when fitting (i.e. models with 
larger cores and higher f3fu can give rather similar profiles 
to those with lower values of both parameters), the pres- 
ence of a central excess will force r core towards lower values 
and hence decrease Pfit, unless an additional component is 
included in the model to account for the central excess. 

Interestingly, we do not find this to be the case in gen- 
eral, for our analysis. For the subset of our groups with 
two-component fits, the median value of j3ju for the single- 
component fits is 0.47 (i.e. just steeper than for the two- 
component fits). Individually, some groups (e.g. NGC533) 
have a steeper profile when the two-component model is 



used, and some (e.g. NGC4761) have a flatter profile. The 
distinction appears to be that the argument of Mohr et al. 
(1999) applies to systems for which the extended component 
dominates over most of the range of the fitted data. In this 
case, the presence of a central component acts to slightly 
modify the extended component fit, by reducing both r core 
and (3 fu. NGC2563 in Fig.| is such an example. However, 
for systems where the central component is more dominant, 
such as NGC3091 in Fig.[| the single component fit is a com- 
promise between a steeper central component, and a flatter 
extended one, so that the result is to increase pfu, relative 
to the extended component. 

Fig.jio] shows the relationship between the f3fu values 
from one and two-component models for the eight systems 
from our sample with well-fitting two-component profiles. 
The solid line splits the graph into two areas. In the up- 
per left area the two-component fit has a steeper profile 
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Figure 10. Relationship between the single-component Pfu 
value and the extended component (3f it value from the two- 
component fit. The groups shown are the eight systems with the 
best fitting two-component models. 



than the single-component fit, in the lower right area the 
reverse is true. As can be seen, the single-component fits 
lead to overestimates and underestimates of /3fit, relative 
to the two-component results, in equal numbers of cases. 
The two dashed lines delineate the region in which the two- 
component fit differs from the single-component fit by less 
than ±50%. As can be seen the two-component models gen- 
erally have (3 fa values for the extended component within 
50% of the single-component fit. The small nominal errors 
on the single-component j3fu values in the figure are mis- 
leading, since they result from calculating errors on a poor 
fit. 

The slope of the L : T relation for our group sample 
is significantly steeper than the cluster relation. This is in 
contrast to the results of MZ98, who find that the L : T 
relation for their sample of nine groups is consistent with 
the cluster relation. However they had too few points to fit 
to the group sample alone, so they added a large cluster 
sample in order to determine the best fit line. If the L : T 
relation is turning over at a temperature of ~ 1 keV, as is 
suggested by FigM then it is to be expected that the line 
fitted through a combined group and cluster sample would 
not differ greatly from the cluster relation. 

Values of /3 3pec derived by MZ98 for their groups lead 
them to conclude that /3 spec ~ 1, whereas we see evidence 
for a drop in f3 spec for low temperature systems (Fig. |^). 
This difference appears to result from two factors. Firstly, 
four of the nine common groups are found in the region 
(T>1 keV) where our groups are generally consistent with 
P 'spec ~ 1. So this only leaves five systems in which MZ98 
could have noted a drop in /3 3pec . Secondly, our values of 
Pspec appear to be typically about 10% lower than those of 
MZ98. For the nine groups in common, we derive a mean 
value of </3 sp ec> = 0.78 compared to <p spe c> = 0.87 for 
MZ98. Since we use the same velocity dispersions, the dif- 



ference results from the derived gas temperatures. This dif- 
ference may arise from the fact that for most groups MZ98 
extract their spectral data from within a larger radius, and 
given the tendency towards a decline in temperature with 
radius apparent in many systems in Figjl], this should result 
in temperatures somewhat lower than ours. This interpreta- 
tion is supported by the fact that our temperatures are in 
good agreement with those derived in the study of Mulchaey 
et al. (1996), in which similar extraction radii were used for 
systems common to the two studies. 



7 DISCUSSION 

This survey of X-ray bright, loose galaxy groups represents 
the largest detailed study of their properties to date. This 
allows a comparison with the properties of richer clusters, 
and we have been able to show that three effects are appar- 
ent in low temperature systems: steepening of the L : T rela- 
tion, steepening of the <r:T relation (i.e. lower f3 sp ec values in 
groups), and flatter surface brightness profiles in groups. We 
find that the contrary results of MZ98 appear to be due to 
the small size of their sample, coupled with their somewhat 
less sophisticated analysis of the surface brightness distribu- 
tions. 

The general nature of these three departures from clus- 
ter trends are in good agreement with the expectations from 
preheating models, in which energetic winds from forming 
galaxies raise the entropy of intergalactic gas and inhibit 
its collapse into the smaller potential wells of galaxy groups 
(Metzler & Evrard 1994; Cavaliere, Menci & Tozzi 1997; 
Cavaliere, Menci & Tozzi 1999; Ponman et al. 1999; Met- 
zler & Evrard 1999; Balogh, Babul & Patton 1999). This 
increase in gas entropy primarily acts to reduce the gas den- 
sity in the central regions of low mass systems, flattening 
their surface brightness profiles and reducing their X-ray lu- 
minosity. The enhanced entropy also leads to some increase 
in gas temperature, resulting in a value of /3 S pec less than 
unity. 

The slope of the L : T relation, L oc T 4 ' 9±as , is flat- 
ter than the index of 8.2 ± 2.7 derived for Hickson groups 
by Ponman et al. (1996), however the error from the HCG 
sample was very large, so the difference in slopes is not sig- 
nificant (1.2cr). The present, much more accurate determi- 
nation of the L : T slope, is in excellent agreement with the 
asymptotic relation L oc T 5 derived in the low temperature 
limit by the semi-analytical models of Cavaliere et al. (1997) 
and Balogh et al. (1999). These two treatments make some- 
what different simplifying assumptions about the physics of 
the heating of the intracluster gas, but both obtain similar 
slopes in the limit of isentropic gas (i.e. where shock heating 
becomes negligible). 

This result has to be quite robust to detailed model 
assumptions, since an approximate result L oc T 4 ' 5 A(T), 
where A(T) is the cooling function, is easily derived by 
combining the scaling relations T oc M/R (from hydro- 
static equilibrium), M oc R 3 (for systems virialising at a 
given epoch), p gas oc y 3 / 2 (for constant entropy gas) and 
L oc p 2 gas A(T)R 3 . For bremsstrahlung, A(T) oc T 1/2 , so that 
one obtains L oc T 5 . In practice, at T ~ 1 keV the cooling 
function is flatter than T 1//2 , due to the increasing contri- 
bution of metal lines at low temperatures, and so the ex- 
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pected relation flattens somewhat towards L oc T 4 ' 5 . The 
good agreement between this isentropic result and our ob- 
servations lends strong support to the result of Ponman et al. 
(f999), that the gas entropy tends towards a constant 'floor' 
value, set by preheating, in low temperature systems. 

Within the above picture, the significant scatter seen in 
our L : T relation is expected to be primarily due to dif- 
ferent star formation and merging histories of the groups. 
It has also been shown (Fabian et al. 1994) that scatter in 
the cluster L : T relation is correlated with the strength of 
the emission associated with a cooling flow. Lower temper- 
ature gas (at a given density) has a shorter cooling time, 
and it is apparent from Fig.hl that many of these groups do 
contain cooling flows. Hence some L : T scatter can also be 
attributed to the presence of cooling flows in the sample. 

Another consequence of the effect of galaxy winds is 
that if winds have injected extra energy into the intragroup 
medium then a greater proportion of the energy of the sys- 
tem should be found in this hot gas. However, this extra 
energy could manifest itself in the form of extra thermal en- 
ergy, or higher gravitational potential energy of the gas. The 
models of Cavaliere et al. (1999) and Balogh et al. (1999), 
and the N-body+hydrodynamical simulations of Metzler & 
Evrard (1999), all indicate that for systems with T > 1 keV, 
the energy is taken up in flattening the gas distribution, with 
very little effect on gas temperature. Unfortunately, the sim- 
ulations of Metzler & Evrard (1999) do not extend to lower 
temperatures, but the models of Cavaliere et al. (1999) and 
Balogh et al. (1999) both predict that at T<0.8 keV, sys- 
tems depart rather suddenly from the cluster M : T relation, 
with T flattening out at a minimum value. This must nec- 
essarily happen, since (in the absence of significant cooling) 
the gas temperature cannot drop below the level to which 
it was preheated, since its density will have increased as it 
settles into the group potential. 

The observed <r:T relation for our groups (Figj^), is 
noisy, but there is a rather clear pattern whereby /3 spec ~ 1 
for T>1 keV, but drops to lower values for cooler systems. 
For example, the median /3 spc . c for our nine groups with 
T < 0.8 keV is 0.44. This behaviour is just what the models 
predict for preheating temperatures ~ 0.5 keV. 

The L : a relation for our group sample is slightly flat- 
ter than the cluster relation as determined by White et al. 
(1997), although the errors on the slope of the loose group 
sample are large and as a result the difference is not sta- 
tistically significant. This might suggest that the group L:a 
relation is an extension of the cluster trend. However, if as 
argued above, preheating has substantially reduced the lu- 
minosity of the groups, then the velocity dispersion must 
also be lower than expected, otherwise a steepening of the 
L : a relation, similar to that seen in L : T, would be ob- 
served. 

Bird et al. (1995) have suggested that velocity disper- 
sion should be reduced for lower mass systems due to the 
effects of dynamical friction, which is more effective in lower 
mass systems due to their lower velocity dispersion. Loss 
of orbital energy will lead to a reduction in orbital veloc- 
ity provided that the potential is less steep than a singu- 
lar isothermal potential in the inner regions. This would be 
the case for either a King-like potential, with a flat core, or 
for potentials of the form introduced by Navarro, Frenk & 
White (1997), which tend to p oc r _1 at small radii. How- 



ever it must be remembered that the velocity dispersions 
of the groups in this sample are drawn from three differ- 
ent sources, and may be based on only a small number of 
group galaxies, so that statistical errors are large. Zabludoff 
& Mulchaey (1998) find that when they add the velocities 
of fainter group galaxies to their redshift samples, the veloc- 
ity dispersions they derive may increase by a factor of 1.5 
or more. This is qualitatively consistent with expectations 
from dynamical friction, since the orbits of more massive 
galaxies should decay more quickly, and hence their velocity 
dispersion would drop below that of fainter group members. 

The results on the asymptotic slope of the X-ray surface 
brightness in groups derived here, confirms and quantifies 
the result of Ponman et al. (1999), who showed that sur- 
face brightness is progressively flattened in low temperature 
systems. This trend is in accord with preheating models, as 
discussed above, although our median value of /3f it = 0.46 
is a little lower than the values /3/jj ~ 0.5-0.6 predicted by 
the models of Metzler & Evrard (1999) and Cavaliere et al. 
(1999) for T ~ 1 keV. 

The situation in clusters is still a matter of debate. Ar- 
naud & Evrard (1999) collect together results from the liter- 
ature, and find a clear trend in /3fu with temperature, as can 
be seen in Fig^. However, Mohr et al. (1999) find that two- 
component fits are required to adequately represent most 
cluster profiles, and that the results from such fits show no 
trend in the value of /3/u for the extended cluster compo- 
nent. They conclude that results such as those of Arnaud 
& Evrard (1999) arise from biases due to the inappropriate 
use of single /3-model profiles. On the other hand, we have 
accounted for the central component, but still find that f3fu 
is substantially lower in groups that the value of 2/3 found 
for clusters by Mohr et al. (1999). 

The resolution of this situation probably lies in the 
temperature ranges covered. The analysis of Ponman et al. 
(1999) is model-independent, in that it involved simply over- 
laying the scaled surface brightness profiles. This shows that 
flattening of the profiles sets in at temperatures T<3 keV. 
Since the sample of Mohr et al. (1999) includes only a single 
cluster with T < 3 keV, the lack of trend in /3/it observed 
within their sample, and the much flatter profiles observed 
in our sample, are both consistent with the Ponman et al. 
(1999) results. 

Finally, we wish to emphasize that an important im- 
plication of the flat X-ray profiles of groups, coupled with 
their generally low surface brightness compared to clusters, 
is that one must be very careful in drawing conclusions about 
properties such as gas mass, gas fraction etc. on the basis of 
analyses confined to 'detection radii'. For example Mulchaey 
et al. (1996) conclude that masses of gas in groups are typ- 
ically lower than the mass in galaxies, on the basis of anal- 
yses within the region of detectable X-ray emission, which 
in many cases is only ~ 200 kpc. Such results have im- 
portant implications. For example, Renzini (1997) has used 
them to argue that the iron mass to light ratio in groups 
is much lower than that in clusters, and that it is therefore 
difficult to explain how clusters can be assembled through 
group mergers. 

It can be seen from Fig|] that under the assump- 
tion that our /3-model fits can be extrapolated to Rv, less 
than 50% of the X-ray luminosity of the system is con- 
tained within 200 kpc for typical groups. Now the asymp- 
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totic power law behaviour of surface brightness at large r 
is S(r) oc r 1-6/3 , whilst the corresponding density profile 
(in the approximation of isothermal gas) is p ga s oc r -3 ' 3 . 
Hence the density profile is even flatter, and the fraction 
of the total gas mass contained within r =200 kpc will be 
considerably less than 50%. The flat gas profiles mean that 
the gas fractions of groups rise strongly with radius, so that 
very different results might be obtained if our instruments 
were sufficiently sensitive to detect group emission out to 
Rv, a possibility which should be realised with the launch 
of XMM. 



8 CONCLUSIONS 

We have carried out detailed analysis of ROSAT PSPC data 
for 24 X-ray bright galaxy groups. Temperatures and bolo- 
metric luminosities have been derived for each group, and 
surface brightness profiles modelled in some detail. In agree- 
ment with previous studies we find evidence for the presence 
of two components in the surface brightness profiles of many 
of the groups. When present, the central component is coin- 
cident with the position of a central galaxy, suggesting that 
it may be due to the halo of the galaxy, or to a cooling flow 
focused onto the central galaxy. 

The surface brightness profiles of groups are signifi- 
cantly flatter than those of galaxy clusters. For a subsample 
of the groups with the best data, the steepness of the sur- 
face brightness profiles, as measured by the parameter 
appear to show a trend with mass when combined with clus- 
ter data. This result is consistent with the idea that galaxy 
winds have significantly affected the state of the intergalac- 
tic medium in low mass systems. 

The relation between the X-ray luminosity and temper- 
ature for galaxy groups is also derived. This relation is found 
to be significantly steeper than that derived for galaxy clus- 
ters. The action of galaxy winds flattening surface bright- 
ness profiles would reduce the luminosity of the gas, due to 
the luminosity dependance on the square of the density, thus 
causing a steepening of the L : T relation for lower mass sys- 
tems. Further evidence for this scenario is provided in the 
relation between velocity dispersion and temperature. The 
a : T relation shows that for lower mass systems the specific 
energy in the gas is greater than the specific energy in the 
galaxies, suggesting that there has been energy injection in 
these systems. An encouraging level of agreement is appar- 
ent between our results and recent models and simulations 
of the effects of preheating by galaxy winds. 
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